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Abstract 



In this paper, it is shown how a single stripe and a stripe phase grow from individual holes in low 
doping regime. In an effective low-energy description of the t-J model, i.e., the phase string model, a 
hole doped into the spin ordered phase will induce a dipolar distortion in the background [Phys. Rev. 
B67, 115103 (2003)]. We analyze the hole-dipole configurations with lowest energy under a dipole- 
dipole interaction and show that these holes tend to arrange themselves into a regular polygon. 
Such a stable polygon configuration will turn into a stripe as the number hole-dipoles becomes 
thermodynamically large and eventually a uniform stripe state can be formed, which constitutes an 
energetically competitive phase at low doping. We also briefly discuss the effect of Zn impurities on 
£*\) ■ individual hole-dipoles and stripes. 



PACS numbers: 71.27+a,71.10-w 



I. INTRODUCTION 



The stripe phenomenon is one of many interesting and novel properties observed in high-T c cuprate superconductors. 
Static stripes were first experimentally found in Lai.4sSro.i2Ndo.4Cu0^i by neutron scattering, where narrow elastic 
magnetic superlattice peaks located at (7t(l±2a;),7r,0) and charge-order peaks at (47r(l±a;),0, 0) are clearly identified 
at doping concentration x — 0.118. This result is interpreted as that the dopant-induced holes collect in domain 
walls that separate antiferromagnetic (AF) antiphase domains. This picture is also supported by the x-ray diffraction 
experiments 2 . In L^-zSr^CuO^ compounds, people have also tried to use dynamical stripes to explain the observation 
by the inelastic neutron scattering^, where narrow magnetic peaks at the AF wave vectors, (tt(1 ± e),7T, 0) and 
(tt, 7t(l ± e),0) with e ~ 2x at low energies were found. Somewhat similar incommensurate dynamic magnetic 
fluctuations were also reported^ in YBCO compounds. Nuclear quadruple resonance(NQR)^, muon spin resonance^, 
and magnetic susceptibility measurements 7 - all verify the evidence of stripe in L^-rrSr^CuO^. Stripes have also been 
observed in the oxygen-doped La2Cu04 using nuclear magnetic resonance (NMR) techniques^. These experimental 
■^- , results suggest that the stripe instability may be extensively present in the cuprates as a competing order, which 
VO contributes to the complexity of the phase diagram. 

The existence of the stripes in a strongly correlated electron system was actually first predicted^ by Zaanen and 
Gunnarsson before the experimental discovery. They found the stripe mean-field solution in a two-band Hubbard 
model, in which holes doped into the parent antiferromagnet generally tend to arrange themselves into straight lines 
aligned parallel to each other, i.e, the charged stripes. Since the experimental discovery of the static stripes in the 
cuprates, theoretical investigations of the stripe and stripe-related physics have been conducted very intensively in 
the high-T c field. Numerical studies of the t-J model by the DMRG present conflicting conclusions as to the existence 
of stripe phases in its ground state, which might be caused by the strong finite-size effects^. Recent theoretical 
developments in the stripe physics have been reviewed in Refsiii. 

So far the most of theoretical studies on the origin of the stripes in the cuprates are either based on phenonrenological 
theories or focused on the static ones at the mean-field level. To truly understand the microscopic origin of the stripe 
phase and its competitive relation with the homogeneous phases (including superconductivity), one needs to know 
when a stripe can be melt and be broken into pieces, namely, what it is made of, and when it can become stable 
against various kinds of fluctuational effects 



Recently, it has been showni&iS that there exists a more stable elementary object, known as a hole dipole, in 
the low-doping spin ordered phase described by the t-J model. Such a charge +e entity can be regarded as a dipole 
composed of a charged vortex (centered at a spinless holon) and a neutral antivortex which is self-trapped in real space. 
Due to the so-called phase string effect, an infinite (logarithmic divergent) energy is needed if one is to "destroy" such 
a composite by separating two poles of the dipole infinitely far away. 

On the other hand, since each hole dipole is self-trapped in real space, its kinetic energy is suppressed. Thus 
the potential energy (from impurities, for instance) and the dipole-dipole interaction between two holes will become 
dominant. In the absence of disorder or impurities and without considering the long-range Coulomb repulsion, 
an inhomogeneous instability has been found in such a system and in particular various stripe instabilities were 
suggested^ to occur. In other words, if a stripe does form in this system, hole-dipoles described above will become 



the elementary building blocks. Consequently the fluctuations and dynamics of stripes as well as the melting of them 
may be understood and mathematically described from a new angle based on the hole-dipoles. 

In this paper, we shall follow up the stripe instability pointed out in Ref^ and demonstrate mathematically how 
a stable stripe can grow from individual hole-dipoles by starting with only a few of them. We find that these finite 
number of hole-dipoles generally form a regular polygon with a minimized potential energy, which is stable against 
the perturbations. With the increase of the hole number, the polygon eventually evolves into a stripe and then stripes 
as the hole concentration becomes finite in the thermodynamic limit, which result in the stripe phase. We further 
consider the Zn impurity effects on both hole dipoles and stripes and predict that stripes can be easily destroyed in 
the presence of random zincs. 

II. THE MODEL 

A. The phase-string model 

We start with the two-dimensional t-J model. At half-filling, it turns into the Heisenberg model with a good 
description of the magnetic properties in the cuprates. The Marshall sign rule is foundi& in the ground state of such 
a model, where the flips of two antiparallel spins at opposite sublattice sites are always accompanied by a sign change 
in the ground-state wave function: ti - * (~ 1) IT- Upon doping, however, this Marshall sign rule will get frustrated 
by the motion of the doped holes. When a hole hops from site to site, a sequence of ± signs will be left behind which 
cannot be repaired by spin- flip processes and is called phase-string 15 . 

The phase-string theory is developed to accurately handle the phase-string effect in the t-J model. This theory is 
based on a new kind of slave-particle formula in which electron operator readsi^ 

Ci a = htb ia (-aye i§i ', (1) 

>era1 
following no double occupancy constraint 



where h\ is a bosonic "holon" creation operator and bi a is a bosonic "spinon" annihilation operator, satisfying the 



hthi + ^2blbia = l. (2) 

a 

Here the nonlocal phase factor e ia precisely keeps the track of the singular part of the phase string effect as well as 
the fermionic statistics of the electron operator, as defined by 

e i^ = e m-^], ( 3 ) 

with 
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and 


^ = ^Imln(^-^)^. 


The effective phase-string model of the t-J Hamiltonian is given by 


H eff =-t h J2[(e^)hlh J+ H. c :_ 


-*£[(. 



(4) 



(5) 



(6) 



(ij) (v)° 



with th ~ t, J s ~ J. The most important and unique structure of the phase-string theory is a mutual dual relation in 
(p)|): For holons, a spinon simply behaves like a ±7r flux-tube and for spinons, a holon also behaves like a n flux-tube, 
which are described by the lattice gauge fields Af • and A^ as follows: 



a.leC 



and 



E4^5>- (8) 

c lee 

for a closed path C with n h liy and nj 1 denoting spinon and holon number operators, respectively. 



B. Holes as dipoles 

In the phase-string theory, the spin-flip operator is defined as 

S+ = (-l) i 6| T 6 i ;exp[^], (9) 

with 

^ = ^2lmhi(zi -zi)n^ . (10) 

In the AF spin ordered phase, the spinons are Bose condensed, i.e., < bi a >=/= 0, with the spins lying in the xy 
plane. The polarization direction of the spin ordering is determined by (Sf) = (— l) 4 < &L >< bn > cxp [«$'']. 
From this, we can see that besides the sign (— l) 4 , which reflects the staggered AF order, there is an additional phase 
exp [«$''] introduced by holons, which represents a twist of the spins with respect to each holon. Namely, each time 
when one circles around a holon once, a 2tt rotation is found in the direction of the spin ordering. The resulting 
spin configuration is called a meron (spin vortex) [Fig. 2 of Refii 9 -]. A meron costs an energy which is logarithmically 
dependent on the size of the system. And for two holons, the induced spin twists are in the same way such that there 
exists a repulsive interaction between them. In order to remove such an unphysical energy divergence, an antimeron 
should be induced^ near every holon-meron to cancel out the spin twists at large distance. An antimeron is defined 
by: 



bia exp 



(11) 



where -df = lining — z°). Here z° denotes the coordinate of the center of an antimeron labelled by k. As a result 

(S?) - (-1)' (&1 T ) (^)exp [t«* - «*] . (12) 

Define 

0* = **-0*=Imln *^§ (13) 

z x + z k /2 

to describe the spin twist, with Zk = e% + ie v k . Here the meron and anti-meron are centered at ±%, respectively. At 
| r^l >> |ejfe|, one obtains a dipolar twist 

k (zxe k )-r t 

% - |r . |2 W 

The energy cost of such a dipole configuration is given as the following 19 

|e fc | 2 
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|r-e fc /2| 2 |r + e fc /2| 2 



<t = nJ s p s c a', (15) 

u 

which is finite. In the above, p s c a 2 = (b + ) 2 (a is the lattice constant). This meron-antimeron spin configuration is 
called a hole-dipol o 18 ' 19 . The displacement connecting the centers of a meron and an antimeron can be defined as the 
dipole moment here. The dipole moment is determined by the two dimensional Coulomb gas theory by a standard 
KT rcnormalization group (RG) methocUSiiS. Near half filling x — ► we estimate the centers of a meron and an 
antimeron as the lattice constant (|efc|) = ro ~ a. 

At the end of this part, we note that in somewhat different contexts, the concept of hole dipoles has been also 
suggested by different authors through different approaches. It is raised in reference s 20 ! 21 that the doped holes 
introduce a local ferromagnetic exchange coupling between their neighboring Cu , which brings frustration to the 
background antiferromagnetism. The frustrating bond acts like a magnetic dipole. From this picture, they studied 
the suppression of antiferromagnetic correlations by the hole-dipoles, and the magnetic phase diagram was obtained. 
It is also pointed out in reference 2 ^ that the doped holes in the form of hole-spin-polaron interact with each other 






FIG. 1: The variables in formula(16) is shown. For two dipoles indexed by k and k , r kk i denotes the distance between their 
centers, e k and ey denote their moments, and ip k and tpy denote the angles between their moments and the line which connects 
their centers. 

through dipolar potential. It is just the excitation of these spin-polarons that forms the strip o 23 i 24 . In reference^, an 
analogy was given between the doped holes in the AF background with the He impurity in the liquid He . As a result, 
it is also found that the mobile-hole creates a long-range dipolar spin-backflow. And in reference 2 ^, the interaction 
between the holes is more carefully considered, including the long-rang Coulomb interaction, the dipolar potential, 
and the short-range attraction. The competing of these interactions leads to the complicated phase diagram, which 
includes the diagonal stripe. 

III. OPTIMAL CONFIGURATIONS FOR MULTI-HOLE-DIPOLES 

A holc-dipole is self-localized^ in space with the suppression of its kinetic energy. If there is an impurity, such a 
hole-dipole can be easily trapped around, whose effect will be discussed in Sec. IV. In this section, we shall consider the 
impurity-free case, in which the hole dipole can be located anywhere in real space due to the translational symmetry. 
For multi-number hole-dipoles, the dipole-dipole interaction will determine the spatial configuration of them. In the 
following we start with the case of two hole-dipoles first. 

A. Energy-minimal configuration for two hole-dipoles 

For two dipoles well separated from each other, the interaction energy is given byiS 



-trd — d 



2q 2 



' ~ KA 2 



e k ■ e k > - 2 — 



l r fefc'| 2 



2g 2 |e fc ||e fc ,| 

= ; & COS((p k + ip k ,) 

r 2 
~ -2 g 2 ^^cos(^ + ^), (16) 

|rkk'| 

in which the size of the dipoles is fixed by r . The alignments of two dipoles are shown in Fig.l. 

Under the interaction l|l(jf) . two dipoles will adjust their dipole moment directions to arrive at the lowest energy. It 
is easy to see that the condition to minimize their potential energy at a fixed distance is 

ip k + ip k > = 0, modulo 2n, (17) 




FIG. 2: Configuration for three dipoles marked by 1, 2, and 3. The arrows indicate the moments of the dipoles. 



which results in an attractive potential energy, 



vf- 



2q 2 r 2 
Irfcfc'l 2 



(18) 



Two dipoles will then move closer and closer until reach a least distance 2ro = 2 |efe'|, determined essentially by the 
size of the dipoles. When the distance between two dipoles are near 2r , the potential described by l(Tp|l usually is 
no longer correct. However, we shall use the formula (|16|l approximately at \rkk' | ^ 2ro and take the positions of the 
dipoles as continuous variables by ignoring the discrete lattice sites in the following considerations of n-dipole case. 
Then the problem is reduced to a mathematical one to search for the optimum configuration for the n-dipoles which 
interact with each other through |_?61) . under the constraint that the distance between any two dipoles should be no less 
than twice of the average dipole moment |2e/y| = 2ro- 



B. Three and four hole-dipoles 



Now let us consider three dipoles case. In Fig. 2, three hole-dipoles are marked by 1, 2, and 3, respectively. 
Suppose the first and second dipoles are located at (±/i, 0), with the angles between their dipole moments and the 
line connecting them to be ±ip. Then we want to find out what is the optimal location for the third dipole. From 
Fig. 2, it is easy to see that for the dipoles marked by 1 and 3, one has 



tp-O x = V 3 
and for 2 and 3, one has 

<p-6 2 = e 1 + 9 2 

Then one finds 

<p = 02 + di 
Defining the coordinate of the third dipole by (x,y), we get 

v 



<p = arctan( 



x + h 



) + arctan(- 



■X + h 



(19) 
(20) 
(21) 

(22) 



and thus 



tan tp 
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arctan(- 

-2hy 

- y 2 — h 2 



arctan(- 



(23) 
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FIG. 3: For two dipoles 1 and 2 with their centers fixed at (±/i, 0) and the angles between their moments and the line connecting 
their centers to be <p, the track of the center of dipole 3 to optimize the direction of its moment is a circle which passes the 
former two dipoles, centering at (0, —h-ctg<p), with the radius h-cscp. The directions of the three dipoles are along the tangent 
of the circle 
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FIG. 4: When the dipole 4 is added, it should make an optimum configuration with any two dipoles. As a result, it locates on 
the same circle mentioned in Fig. 3. 



such that 



(y 



tamp' 



sin ip 



(24) 



^From (|24(l . we can see that the track of (x,y) is just a circle passing through dipole 1 and 2, centered at (0, —h/ tamp) 
with a radius R = hj sin p. From the knowledge of geometry, it is easy to see that the moments of these three dipoles 
will all point along the tangents of the circle, as shown in Fig. 3. 

Then the energy-minimal configuration for three hole-dipoles is obtained as an equilateral triangle. The minimal 



interaction energy for three hole-dipoles is 



2 



d-d o o„2 r 



V 2 a - a ~ -3 • 2q , 

= -2<Z 2 % (25) 

with |rjwfe/ 1 = V3^ and i? is radius of the circle crossing three dipoles. 

When the fourth hole-dipole is added, it should first make up an optimal configuration with the dipole f and 2, and 
thus is located on a circle as discussed above. Then by further making up an optimal configuration with the dipole 1 
and 3, it will also be located on an another circle passing through the dipole 1 and 3. Since the dipole 3 is already on 
a circle determined by f and 2, with its dipole moment along the tangent, the circle determined by the dipole 1 and 
3 is the same as that determined by the dipole f and 2, and so is the one determined by dipole 2 and 3. 

The same argument is applicable to n dipoles. So the optimal configuration for n dipoles under the potential l|16|l 
will be always on a circle, with each dipole moment along the tangent of the circle and the proportional spacing for 
each dipole. Such configuration is just a regular polygon with n edges. To low energy the radius of the circle R will 
shrink until the length of every edge is equal to the minimum 2ro, 

n 

The minimal interaction energy for n hole-dipoles (n is a odd number) is 

1 



Vt d - -VrgJ] 



kk' 



I 2 

rfcfc' 



2q\t-n-2-Y / -±-2 1 (27) 



with 



|r fc |=2flsin(-fc) = 2r Sm( V . (28) 

n sin - 



,7T 



If n is even number, the minimal interaction energy is 



-l 



V^-^np.g^-L + y^l. (29) 

Numerical simulations for up to n = 30 dipoles also show that for these interacting dipoles, the optimized configu- 
rations with minimized total energy are always regular polygons as discussed above. 

C. Stability of the regular polygon configuration 

In this part we will show the stability of the regular polygon configuration for hole-dipoles. For a regular polygon 
configuration, the distance between two nearest dipoles has reached the minimal 2ro. But the distances between other 
pairs of dipoles have not arrived at their minimum. Mathematically, it should be proved that the configuration is 
stable against perturbations. 

1. Perturbation for changing the direction of a dipole moment 

First we change the direction of a dipole moment from it's original direction ipo to ipo + Sip. The change of the total 
energy for n hole-dipoles (n is a odd number) is 

Er 2 
—^ [cos (ipo + 8ip + ip k ,)- cos (tp + ip k , )] 

k Kl 

A)(M 2 , (30) 

kfcl 
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FIG. 5: When more dipoles are added, they should all locate on the circle mentioned in Fig. 3 and Fig. 4. Here n = 6 is displayed. 



2r-- 



2r„ 



2r„ 



2r- 



2r„ 



2r„ 



FIG. 6: The optimum configuration for n dipoles. Their centers form a regular polygon. The moment of each dipole is directed 
along the bisector of each external angle. Here n — 6 is displayed. 



with |rjfe| defined by 1)28)1 . On the other hand, if n is even number, the energy difference is 

1 1 1 



AV d ~ d ~ q 2 r 2 [2 ■ £ * + ^j^}(Srf > 0. 

fc =i Kl 



2 (2i?,) 2 



(31) 



When more hole-dipoles change the directions of their moments simultaneously, the energy cost is simply the sum 
of all the positive energy cost by every change. 

So the regular polygon configuration is stable against changing the directions of the dipole moments AV d ~ d > 0. 
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FIG. 7: A perturbation is imposed to the optimum configuration, .i.e, the regular polygon. Here n — 4 is displayed. Each 
dipole has a perturbation in the position and direction (Sxi,5yi,89i). The dipole and dipole 1 are fixed at the origin point 
and on the x axis respectively, as the total potential has global SO(2) invariance. 

2. Perturbation for changing the positions of the centers of the dipoles 

Next we change the positions of the centers of the dipoles. In Fig. 7, the case of n = 4 is shown as an example. The 
position and direction of each dipole are marked in Fig. 7. Here 2ro is set to be 1 for convenience. The dipole and 
dipole 1 are fixed at the origin point and on the x axis respectively. Each dipole has a perturbation in the position 
and direction {6xi,6yi,89i). After this perturbation, we can expand the potential up to first order. Taking dipole 1 
and dipole 2 as an example, we have 



among which 



V 12 



-Vocos(ipi + ipi) 



(32) 



'12 



For the numerator, we have 



V = 2q 2 r 2 . 



Vi + ¥>2 = + 8$i + 89 2 - 25012. 



(33) 



(34) 



In the above equation, 89\ and 602 denotes the changes in the directions of the moments of dipole 1 and dipole 2, 
while S9i2 denotes the change in the direction of line connecting dipole 1 and dipole 2, 

<5#i2 = S(x\ — X2)/1 = 5x\ — 8x2- (35) 

Thus the change in the numerator is zero to the first order of the perturbation. For the denominator, we have 



,.2 
12 



/■T , - {6xi - 8x 2 ) 2 + (1 + 5y 2 f 

1 + 25y2 + o{8x 1 , Sx 2 , 6x16x2)- 



As a result, up to the first order perturbation, the potential becomes: 

-V 



V v . 



1 + 25y 2 



= -V + 2V 8y 2 + o{8yl). 



(36) 



(37) 
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In the same way, we expanded the potential among other dipoles, and obtained the following expansion of the total 
potential up to the first order perturbation , 

6Vt d oT d = KVoi + V 12 +V 23 + Vo2 + Vo3 + V u ) 

= -V {6x 1 +5x 2 -5x 3 + 5y 2 +5y 3 ). (38) 

To satisfy the constraint that the distances between the centers of any two dipoles should be no less than 2ro, we 
expanded the formula of their distances to the first order perturbation and have: 

Sxi > — * roi > 1, 
5x 2 > Sx 3 -* r 2 3 > 1, 
Sy 2 > -> r 12 > 1, 

Sy 3 > -► r 03 > 1. (39) 



^From (|38|l and l|39|l . we obtained that: 



SVtoi > 0, (40) 



which denotes that the change of the total potential is no less than zero up to the first order perturbation. 
Our above demonstration can be easily generalized to arbitrary n. For a general n, we have 

V tot =Y. ~ V ° C0S{ ^ +V3) - (41) 

ij y 

Up to the first order perturbation, we have: 

sv = v V 2(yXl - Eil^! I 5x ^ + 2 fa ~ j/j)(% ~ %) ( 42 ) 

The constraint for the least distance between any two adjacent dipoles reads: 

Sri,i+i = 2(xi - x l+ i)(Sxi - 6Xi+i) + 2(y 4 - y i+1 )(Syi - 8y i+1 ) > 0. (43) 

It can be checked that 

6Vtot = gV Y, Sr i,i+i, (44) 

i 

in which 

= J27=2 sm[7r(i - l)/n]/[sin(7ri)/n] 3 
9 8sm(ir/n)sin(2ir/n) [ j 

^From (gSJ, 121) and igSJl, we have 

SVtot > 0- (46) 

^From Ij46(l , we know that the change of the total potential is no less than zero up to the first order perturbation for 
arbitrary n. When all the Sxt and <5y, are carefully chosen so that the "=" is realized in (|43|) . I|46|l turns into 

SV tot - 0. (47) 

In such cases, we should have to check the second order perturbation in the potential energy as the first order 
perturbation is zero. We again gave our proof for n = 4 as an example. 

When 8x\ = Sy 2 = 8y 3 = , and 5x 2 — 8x 3 — 8x, the "=" is realized in l|39l) and hence (|40l ) . The first order 
perturbation of the total potential is zero, so we expanded it to the second order perturbation, and obtained 



Vtoi > E- 

-2Fo 



-Vo 

r 2 - 

1, j n 

21/ V Vo 



l + (5x) 2 l + (l-5x) 2 l + (l + dx) 2 



-hV a + ^V a {5xf+o{{5xf), (48) 
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and therefore, 

SV to i > ~Vo(6x) 2 > 0. (49) 

So the change in the total potential is also no less than zero up to the second order perturbation. 

In the case of n > 4, we can carry out similar expansion, and draw the same conclusion. 

By the above perturbative expansion, we proved that the configuration shown in Fig. 6 is stable against local 
perturbation in the positions and directions of the dipole moments. 

D. Energetically minimized configurations for infinite number hole-dipoles - stripes 

In the above part, we have proved that n dipoles will arrange themselves to form a regular polygon to minimize 
the energy. And the radius of the circle R will shrink until the length of every edge is equal to the minimum Itq. 
When n — > oo, the radius of the circle R turn to diverge R ~ ro IL ~ ^-n — ► oo and a regular polygon will naturally 
be stretched into a line (i.e., a stripe), as shown in Fig. 8. The minimal interaction energy for a regular polygon with 



large number of hole-dipoles is about 

V d ~ d ~ _ 2g 2 r 2 ^ — 



2 
kk> \ rkk '\ 
2„ 



~ -<fn (50) 

When a dipole leaves away from the line ( a stripe), finite energy AV d ~ d will be cost 

AV d ~ d ~ -2q 2 rlY -J-j + 2q 2 r 2 y -L 

= ^(*r) 2 , (51) 

where rj, = 2kro and |r/-/| = \2vk\ + (Sr) 2 . Thus the line shape configuration - stripe is stable against local 
perturbation. 

It is easy to see that such a "stripe" of charge carriers is embedded in a domain wall of the AF background. To see 
this, we consider a stripe along the s-axis composed of the hole-dipoles of a size |e| = r$ and spaced by I = a |e|. Far 
away from the x-axis, the total spin twist summed from i|14|) is given by 



E 



k k ik 

* -sgn(2/i), (52) 

a 

when \yi\ ^$> ro ~ a. Here y^ = yi and Xiu — Ik + x% according to the definition. Thus a phase shift is found across 
the stripe with 

7T 

A(j> = (j)y >0 - (f>y <Q = 2 • — . (53) 

For a special case, a = 2, the line becomes an antiphase domain wall 

7T 

A<p = 4> y> u - 4> y< u = 2 • — = 7r. (54) 

Namely, a stripe composed of hole-dipoles is topologically an antiphase domain wall. 

^From H54p . we can also understand physically the reason why the holes tend to arrange themselves into a straight 
line. When the hole-dipoles are distributed randomly in the AF background, each dipole induces a spin twist as 
described by ((Fiji , which costs additional energy. But when the hole-dipoles form a stripe, the total twist spins away 
from the domain wall will be cancelled out such that the spin ordering on either side of the domain wall becomes 
unfrustrated, just like at half-filling. 

Yet there is a further advantage in the formation of the stripe, i.e, the kinetic energy of the holes can be gained. 
Recall that an isolated hole-dipole cannot move freely as it is self-trapped in real space. But when the hole-dipoles 
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FIG. 8: When n —* oo, the regular polygon shown in Fig. 6 turns into a line. The distance between the centers of neighboring 
dipoles is 2ro. 



2r„ 



FIG. 9: When more holes are doped to occupy a certain concentration of the lattice grid, parallel lines are formed to make the 
stripe phase 

arrange themselves into a straight line, the individual holes actually may move freely along the stripe such that a 
delocalization energy can be gained. This will correspond to a metallic stripe case. 

With a further increase of doping, more stripes will be formed as shown in Fig. 9. For any two adjacent stripes, 
without the spin frustration inside the domain between them, they will not gain additional energy by being closer. 
Under a long-range Coulomb repulsion, a uniform stripe phase will be stable against a cluster formation. In this 
uniform stripe phase, the distances between two neighboring stripes are the same, which obviously is determined by 
the hole concentration. It should be noticed that the recently found "checkerboard" pattern in the LDOS by the 
STM in the cuprate may display a new kind of CDW order—. Its possibility in the framework of the above mentioned 
dipole picture is under further exploration when the long-range Coulomb repulsion is considered more carefully. 



13 



p. o Zn 



FIG. 10: When a hole and a Zn 2+ are doped into the AF background simultaneously, they would like to attract each other 
and form a Zn + -holon dipole. H denotes the holon, while Zn denotes the Zn + 

IV. THE EFFECTS OF ZN-IMPURITIES 

As emphasized before, the stripe formation is intimately related to the self-localization of individual hole-dipoles. 
We also have mentioned that disorder or impurities in the system may have an "amplified" effect on localization. 

In the following we discuss the effect of Zn impurities on hole-dipoles as well as the stripe phase. It is well known 
that doped Zn atoms will be present in the form of Zn with a closed-shell structure, substituting the Cu sites in 
the Cu02 planes of the cuprates. In the t-J model, the site occupied by the Zn may be imposed by a boundary 
condition of an "empty" site where no electron or hole can stay there at low energy. 

Let us examine how a Zn impurity and a hole-dipole will interact. Inside a hole-dipole, those spins on a loop circling 
the center of the antimeron will have a 2ir rotation in their polarization directions. As the radius of such a loop shrinks 
continuously to the antimeron core, the spin polarization directions will change quickly and becomes uncertain at the 
core site. So a spin at the core site of the antimeron will just like a "defect" spini^, and the bonds which connect 
the core spin with its surrounding spins can be thus viewed to be effectively "cut off', resulting an energy increase 
roughly ~ 4 J' (J 7 denotes the average superexchange energy for one bond). On the other hand, when a Zn 2+ is doped 
and replaces a normal Cu site, the bonds which used to connect such a Cu with its surroundings are cut off with 
an energy cost approximately ~ 4J'. But, if such a Zn is to replace the Cu sitting at the core of the antimeron 
of a hole-dipole, no additional superexchange energy will be cost by breaking up those four bonds connected to the 
impurity site. Therefore, it will be energetically favorable for a hole-dipole to be trapped by a Zn impurity. As shown 
in Fig. 10, a Zn impurity locates at the center of the antimeron. Experimentally, there is evidence^ that doped holes 
are indeed trapped by the Zn impurities. 

According to the above discussions, a large number of Zn impurities will not favor the stripe formation, for they 
tend to trap hole-dipoles around themselves. The random distribution of the Zn leads to the random distribution 
of the hole-dipoles. If the Zn concentration is very low, then a stripe is expected to be easily pinned by a Zn, and 
be bent in order to pass several zincs, as shown in Fig. 11. From the transport measurements^, and also from the 
muon-spin- relaxation (|ISR) measurements^, it is found that a small amount of Zn impurities are effective for the 
pinning. So we predict that Zn impurities are very effective in destroying a stripe phase at low doping. 

V. CONCLUSION 

In this paper, it is shown that in the framework of the phase-string model, each hole doped into a spin ordered 
phase at low doping will act as a dipole which is self-trapped in real space. With the suppression of the kinetic 
energy, the dipole-dipole interaction between hole-dipoles will dominate the low-energy physics in the absence of 
disorder or impurities. We demonstrated in detail how a few hole-dipoles collapse into a stable configuration of the 
regular polygon, which turns into a stripe (stripes) in thermodynamic limit. Consequently, we found that at a finite 



14 



Zn Q ..._...._^.._^ e Zn 



\ 



\ 



.0 
Zn 



FIG. 11: A stripe in the presence of the Zn 2+ . The stripe has to pass the Zn 2+ , and so is pinned. 

concentration of holes at low doping, a uniform stripe phase is highly competitive. The effects of impurities are also 
discussed. When a Zn + is doped into the system , it generally tends to trap a hole around itself to form a Zn + -holon 
dipole. As a result, the stripe will be pinned near the Zn site. Further more, we predicted that a finite concentration 
of zincs can easily kill a uniform stripe phase. 

Acknowledgments 

We acknowledge stimulating discussions with Z. Y. Weng. This work is partially supported by NSFC grant no. 
90103021, no. 10247002, and no. 10204004. S. P. Kou acknowledges the support from Beijing Normal University. 



2 



Electronic address: fyang@castu.tsinghua.edu.cn 

J.M.Tranquada, et al, Nature.375, 561(1995); J.M.Tranquada, et al, Phys. Rev. B 54,7489(1996). 
Bianconi, Phys. Rev. B 54, 12018; M.V. Zimmermann, et al, Europhys. Lett. 41, 629 (1998) 

3 S-W. Cheong, et al, Phys. Rev. Lett. 67, 1791(1991); T.E.Mason, et al, Phys. Rev. Lett. 68, 1414(1992); K.Yamada, et al, 
Phys. Rev. Lett. 75, 1626(1995). 

4 P. Dai, et al, Phys. Rev. Lett. 80,1738(1998). 

5 F.C. Chou, et al, Phys. Rev. Lett. 71, 2323(1993). 

6 F. Borsa, et al, Phys. Rev. B 52, 7334 (1995). 

7 J.H. Cho, Phys. Rev. Lett. 70, 222(1993). 

8 B. W. Statt, Phys. Rev. B 52, 15575 (1995); N. J. Curro, et al, Phys. Rev. Lett. 85, 642(2000). 

9 J.Zaanen, O.Gunnarsson, Phys. Rev. B 40, 7391(1989). 

S.R. White and D.J. Scalapino, Phys. Rev. Lett. 80, 1272 (1998); 81, 3227 (1998); 84, 3201 (2000); C.S. Hellberg and E. 
Manousakis, Phys. Rev. Lett. 83, 132 (1999); 84, 3022 (2000). 

1 J.Zaanen, Nature 403, 714(2000); Science 286, 251(1999) 

2 G. Kotliar and J. Liu, Phys. Rev. B. 38, 5142 (1988) 

3 D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 (1988) 

4 Shraiman and E. D. Siggia, Phys. Rev. Lett. 62, 1564 (1989); 61, 467 (1988), C. Jayaprakash et al, Phys. Rev. B 40, 2610 
(1989); D. Yoshioka, J. Phys. Soc. Jpn. 58, 1516 (1989); C. L. Kane et al., Phys. Rev. B 41, 2653 (1990). 

5 Z. Y. Weng, et al, Phys. Rev. B. 55, 3894(1997) 

6 W. Marshall, Proc. R. Soc. London Ser. A 232, 48 (1955) 

7 Z. Y. Weng, et al, Phys. Rev. B. 59, 8943(1999) 

8 Su-Peng Kou and Zheng- Yu Weng, Phys. Rev. Lett. 90, 157003(2003) 



15 



1H S. P. Kou and Z.Y.Weng, Phys. Rev. B 67, 115103 (2003). 

20 A. Aharony, et al, Ph ys. Rev. Lett. 6 0, 1330(1988). 

21 A. Cherepanov, et al, |arXivx ond-mat/9808235| 

22 J.R. Schrieffer, X.G. Wen, and S.C. Zhang, Phys. Rev. B 39, 11663(1989); D.M. Frenkel and W. Hanke, Phys. Rev. B 42, 
6711 (1990). 

23 A. Singh, Z. Tesanovic, Phys. Rev. B. 41, 614(1990); 41,11457(1990) 

24 A. L. Chernyshev, et al, Phys. Rev. B. 65, 214527(2002) ;Phys. Rev. Lett. 84, 4922(2000) 

25 Efstratios Manousakis ,Nucl. Phys. B 20, 689(1991) 

26 Branko P. Stojkovic, Phys. Rev. Lett. 82, 4679(1999), Phys. Rev. B 62, 4353(2000), arXiv:cond-mat/9911380 

27 J.E. Hoffman, et al, Science 295,466 (2002); Science 297, 1148 (2002); C. Howald et al, Phys. Rev. B 67, 014533 (2003); M. 
Vershinin et al, 10. 1126/science. 1093384 (Science Ex- press Reports), February 12, 2004; H C Fu , J C Davis and D-H Lee, 
cond/mat 0403001; P.W.Anderson, cond/mat 0406038; H Chen et al., cond/mat 0402323. 

28 M. Hucker, Phys. Rev. B. 59, R725(1999) 

29 Y. Koike, et al, J. Low Temp. Phys. 105, 317 (1996); Y. Koike, et al, Physica C 282-287, 1233 (1997); T. Adachi, et al, J. 
Low Temp. Phys. 117, 1151 (1999). 

30 I. Watanabe, and K. Nagamine, Physica B 259-261, 544 (1999); I.Watanabe, et al, J. Phys. Chem. Solids 63, 1093 (2002); 
I.Watanabe, et al, Phys. Rev. B 65, 180516(R) (2002); I. Watanabe, et al, Physica B 326, 305 (2003); T. Adachi, et al, J. 
Low Temp. Phys. 131, 843 (2003); T. Adachi, et al, cond-mat/0306233 



